Self-Organized Ordering of Nanostructures Produced by Ion-Beam Sputtering 
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We study the self-organized ordering of nanostructures produced by ion-beam sputtering (IBS) of 
targets amorphizing under irradiation. By introducing a model akin to models of pattern formation 
in aeolian sand dunes, we extend consistently the current continuum theory of erosion by IBS. We 
obtain new non-linear effects responsible for the in-plane ordering of the structures, whose strength 
correlates with the degree of ordering found in experiments. Our results highlight the importance 
of redeposition and surface viscous flow to this nanopattern formation process. 
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Performance of many of the (opto)electronic devices 
currently being designed based on arrays of nanostruc- 
tures such as quantum dots, requires a high degree of 
in-plane ordering pj . Currently, there is a formidable ef- 
fort to develop experimental techniques which are able to 
provide highly ordered nanostructures in a self- organized 
fashion 2]. These would allow for easy, low-cost and 
large area fabrication of patterned structures. Among 
these techniques, erosion by ion beam sputtering (IBS) 
at low energies Q is a promising candidate Q |, ||, 
leading to production of nanostructures of varying de- 
gree of uniformity and order, onto diverse substrates such 
as GaSb, InP, and Si. Therefore, detailed knowledge of 
the basic mechanisms underlying erosion by IBS is cru- 
cial to understand and control the associated manufac- 
turing process. From a fundamental point of view, the 
dynamics of surfaces eroded by IBS exemplifies neatly 
the interplay of fluctuations, external driving and dy- 
namic instabilities, sharing many features with seemingly 
unrelated non-equilibrium systems, such as aeolian sand 
dunes 0,E|- Thus, typically surfaces eroded by IBS spon- 
taneousiydevelop submicrometric patterns (dots, pits, 
ripples) [9j depending on experimental conditions, may 
deteriorate and eventually lead to rough interfaces, with 
fluctuations described by the universality classes of ki- 
netic roughening [lfj|. 

A successful approach to surface erosion by IBS is pro- 
vided by continuum evolution equations for the surface 
height, allowing access to time and length scales typi- 
cal of the corresponding pattern formation process. This 
approach was pioneered by Bradley and Harper (BH) 
[Llj, who, based on Sigmund's linear cascade approxima- 
tion of sputtering in amorphous or polycrystalline tar- 
gets , derived a linear equation that describes satisfac- 
torily the main features of ripple formation under IBS, 
such as their alignment with the ion beam as a function 
of incidence angle. Additional features, such as ripple 



stabilization, wavelength dependence with ion energy or 
flux, or production of dot or hole structures as a function 
of bombardment conditions, required extensions of BH's 
approach 0, 0] , leading to a non-linear equation of the 
Kuramoto-Sivashinsky (KS) type. The KS equation pro- 
vides the continuum description of interfaces appearing 
in many diverse systems, see in 12], in which a periodic 
pattern develops, with a preferred wavelength (i.e. the 
lateral size of the nanostructures), that evolves into a 
disordered array. Thus, the crucial properties of homo- 
geneity and in-plane short-range hexagonal ordering of 
the nanostructures produced by IBS remain to be under- 
stood • Recent attempts have been made at exten ding 
the KS equation to overcome such shortcomings [lfll llfij . 
that do not provide definitive answers, since they either 
conflict with symmetries of the physical system, or with 
mathematical requirements for well-posedness |l7| . 

In this Letter, we present a new continuum model 
of erosion by IBS. It leads to a physically and math- 
ematically well defined generalization of the KS equa- 
tion that explains within an unified framework the vary- 
ing degree of homogeneity and order of the nanostruc- 
ture arrays produced 0, 0, @> as a function of exper- 
imental parameters. We exploit connections with rip- 
ple formation in sand dunes Q, hinted at by Aste and 
Valbus a 1181 . o vercoming limitations of previous theories 

Model- During IBS, the bombarding ions penetrate 
the target and induce complex collision cascades in the 
bulk. In semiconductor substrates like those studied in 
, these cascades amorphize the near-surface layer. 
Sputtering events take place when surface atoms receive 
enough energy and momentum to break their bonds and 
leave the target. We will assume that only a fraction 
of those atoms are redeposited at the surface. Adatoms 
are moreover available to relaxation mechanisms such as 
surface diffusion, that can be thermally activated, or else 
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be induced by the mentioned change in the local viscosity 
of the material close to the surface . 

In the spirit of the so-called hydrodynamic theory of 
ripples in aeolian sand dunes Q, we define two coupled 
fields, namely, R(x,t) and h(x,t), where x = (x, y). The 
first one represents the fraction of surface atoms that are 
not sputtered away but, rather, remain mobile along the 
target surface. Analogously, h measures the height of the 
surface neglecting the contribution from the fraction of 
mobile atoms R. Time evolutions of R and h are coupled 
through reaction and transport mechanisms H3. Thus, 
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where T ex and T ac i are, respectively, the rates of excava- 
tion and addition to the surface, v is the average velocity 
of mobile atom, and 4> is the fraction of adatoms that 
detach irreversibly from the surface. Thus, system JU-01 
does not conserve the amount of material, in marked con- 
trast with typical conditions for aeolian sand dunes ||. 
Here, large redeposition of sputtered atoms corresponds 
to the small <f> limit, while, in the absence of redeposi- 
tion, = 1. Considering that matter transport along 
the surface is due to diffusion of mobile species, we set 
— V • J = DV 2 R, where D is the surface diffusivity. 

In the absence of bombardment, the concentration of 
mobile adatoms R changes due to thermal nucleation 
of adatoms from the "immobile state" h, and subse- 
quent transport along the surface. Assuming nucleation 
events are more likely in surface protrusions, we have 
r Sd° r ' = T ^[ R - Reqi 1 + Ak )]> analogous of the Gibbs- 
Thompson relation, n being the mean surface curvature 
and A the capillary length, assumed isotropic due to 
amorphization by the ion beam. Here r is related to 
the mean time between nucleation events, and R? eq is the 
mean equilibrium concentration of mobile species for a 
flat surface. In the presence of bombardment, r"^ cr ' has 
to be generalized, to include the contribution of erosion 
to surface mobility [l^. If the ions fall onto the target 
along the x direction, forming angle 9 with the normal to 
the uneroded target, we have, for small slopes [2(1 EH] . 

Y ex = a [l + Ai 2 (V/i) 2 ](l + <*i • Vh + a 2 W 2 h) 

- a [a 3 {Vh) 2 -a 4 (d x h)(V 2 h)]-pV 2 h, (3) 
T ad = 7o [i?-i? e9 (l- 7 2V 2 /i)], (4) 

where R eq and ji generalize parameters in r"^ cr ' so that 
7o = t -1 + T~J-, 72 = A + A ex , with t~ x and A ex be- 
ing analogs of nucleation time and capillary length of 
erosive origin 0, ■ Coefficients a* > in (J3J are re- 
lated to geometric correction factors that take into ac- 
count the local variation of the ion flux with the surface 
slopes E.g., for oblique incidence, a 1,0:4 oc sin#, 

and 03 = 1/2. Likewise, coefficient fi 2 > is related to 
the local variation of the sputtering yield with the surface 



slope |22|, assumed to have a local minimum for normal 
incidence, while (3 > measures the efficiency of erosion 
due to direct impingement of the ions onto surface atoms 
(knock-on sputtering) 0, ■ The positive sign of a 2 im- 
plements the physical instability inherent to Sigmund's 
theory, by which erosion is more efficient at surface de- 
pressions than at surface protrusions Q. Actually, the 
analysis presented below will allow us to relate some of 
these coefficients with the parameters characterizing Sig- 
mund's distribution of energy deposition. 

Surface dynamics.- Our continuum model of IBS, 
formed by Eqs. C3>-© ; ©-©) provides a way to intro- 
duce systematically all relevant physical mechanisms for 
IBS, differing from that in 0] in a number of features. 
Rather than considering its full solution, we proceed by 
deriving an effective equation for the surface height. As 
in the experiments of references 0, 0, @i we consider 
the case of ions bombarding the target at normal inci- 
dence (0 = 0), thus ai = a 4 = |v| = in ©, © 
After a transient time of order Jq 1 , Eqs. have 
a planar solution ho(t) = —ao4>t, Ro(t) = R eq + (1 — 
(j))ao/jo. Perturbing this solution with periodic waves 
of the form hk = h% exp(^fet + ik • x), and an analogous 
expression for Rk, amplification/decay of such pertur- 
bations is characterized by the dispersion relation ojk = 
Re q lol2 (e# 2 - Jo\D + cf>A)[l - e(l - 4>)}k 4 ) , with e = 
A/(i? eg 7o72) and A = aoa 2 — (3. If A > in cjfc, i.e. if 
sputtering is dominated by collision cascades rather than 
knock-on events, as occurs at low to intermediate ener- 
gies where Sigmund's theory is applicable, there is a band 
of unstable modes that grow exponentially fast, with a 
linear dispersion relation cut of the expected KS type. At 
this stage, the surface morphology is dominated by a peri- 
odic pattern whose wave- vector maximizes u>k- In-plane 
isotropy under normal incidence implies dependence of 
LUk on k = |k| rather than the full wave- vector k, thus the 
surface power spectral density is, rather, maximum on a 
ring [a, |l4| . Stabilization of this pattern occurs when its 
amplitude is large enough that non-linear effects are no 
longer negligible. Close to the instability threshold, the 
rate of erosion is much smaller than the rate of addition 
to the surface. Hence, parameter e above, which is the 
ratio between these two typical rates, is small. We thus 
can perform a multiple scale expansion by introducing 
time scales T% — et and T 2 — e 2 t, and by rescaling length 
scales as X — e 1 / 2 ^. To lowest non-linear order 0(e) and 
as seen in the slow variables (X. T = T\ + T 2 ), surface 
dynamics is described by (see |2fJ for details) |24| 

d T H = -isV 2 H - JCV 4 H + Ai(ViT) 2 - A 2 V 2 (ViI)( 2 q) 

where H = h\ + eh 2 , and 

v = A<f>, 1C = e% l {D + <pA)[R eqlQl2 ~ A{l-<j))], 
Ai = 0oo(l/2-Ma), (6) 
A 2 = e ao(l/2- M2 )7o _1 [(^ + ^)(l-0)-^e 9 07o72]- 
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Eq. JSJ with a noise term, has been already employed 
in the growth of amorphous thin films |25|. In our con- 
text, Eq. JSJ has some important limits. First, in the 
absence of ion bombardment, A = a = 0, 70 — » t _1 
and 72 — * A, and in the original variables (JSJ reduces 
to Mullins' equation for thermal surface diffusion [iif . 
9 t ft- = -L»i? e ,AV 4 /i. In the general case, (J5J-|JBJ| include 
contributions to surface diffusion that arc both thermally 
activated, and directly induced by the ion beam as in 
|lf)j . Second, the BH limit corresponds to 4> — 1, i.e., 
no redeposition. While in |18| the BH limit zeroes out 
the k A contribution to the analog of ujk — thus making 
the typical length scale of the dot structures remain un- 
defined within linear instability — , here Eq. J^} recovers 
for <f> = 1 the equation obtained within BH's approach to 
Sigmund's theory [l2L Hr3 ] , including the fact that the co- 
efficients of the two nonlinear terms have the same signs 
thus making the equation nonlinearly unstable and math- 
ematically ill-posed 0,0] . Thus, beyond its experimen- 
tal relevance, redeposition is crucial in order to make the 
theory mathematically sound. On the other hand, the 
BH limit allows us to extract the phenomenological de- 
pendence of the parameters in our model with character- 
istics of the collision cascades, such as the ion penetration 
depth, a, and the longitudinal and lateral widths a, and 
fi, characterizing the Gaussian decay of enery deposition 
|!|. Thus, for 4> —> 1 we have, in the notation of [l^ . 
a = F,a 2 = a^ 2 /(2a 2 ), ^ = 1 - fi 2 /(2a 2 ) - ^ 2 /{2a A ), 
Reqj2 = A* 2 /4, with F oc JE/cr, where J and E are the 
average ion flux and energy, respectively. 

Eq. (JjJJ describes the evolution of the erosion process. 
Initially, dynamics is controlled by the linear terms, with 
the same dispersion relation ujk as above, and a periodic 
pattern develops, with characteristic wavelength given by 

l c = 2tt [2R eql2 {D + <i>A)[l - e(l - 0)]/(A/>)] 1/2 , (7) 

providing the typical size of the nanostructures that form. 
When local slopes become large, the nonlinear terms in 
Eq. JSJ control the dynamics in an opposing way. While 
the A2 term tends to coarsen the nanostructures in ampli- 
tude and lateral size, similarly to its role in the coarsening 
of ripples on aeolian sand dunes Q, the nonlinearity Ai 
tends to disorder the pattern leading to the paradigmatic 
KS spatiotemporal chaos. Remarkably, Ai(V/i) 2 seems to 
interrupt the coarsening process induced by — A2V 2 (V/i) 2 
and the stationary state morphology consists of domains 
of hexagonally ordered nanostructures separated by de- 
fects. The density of these is a function of the ratio 
r = A2/A1, whose r — > limit in Eq. (JSJ leaves us with the 
KS equation. In Pig. W[a), we plot the stationary-state 
morphology obtained by numerical integration of Eq. (JSJ 
for a relatively large ratio r = 5 27] . The high degree of 
in-plane short range hexagonal ordering is made clear by 
the height autocorrelation function, shown in the inset 
of Fig. Ofa) . The time evolution of the dot pattern can 



be assessed in Fig. |2Ja) , in which the surface roughness 
(mean height square deviation) W{t) vs t is shown for 
the same parameters as in Fig. [T^ al . In excellent agree- 
ment with measurements for nanodots on GaSb , the 
roughness first increases exponentially during develop- 
ment of the linear instability, attains a maximum value 
after dots have coarsened to form a densely packed array, 
and finally relaxes to a smaller stationary value when de- 
fects among different dot domains arc annihilated. Times 
between linear instability and maximum in the roughness 
correspond to non-linear coarsening of the dot structures, 
as seen in the plot of the lateral correlation length £ c (i), 
shown on the same panel. We define £, c (t) as the length- 
scale provided by the first secondary maximum of the 
height autocorrelation. As seen in Fig.Efo), £ c (i)scale is 
constant during linear instability, grows as i°- 27±0 02 j and 
saturates at long times, in agreement with experiments 
on InP 5] . This interrupted coarsening process has been 
also observed on Si 6] (6) and GaSb Q. 

Experimental conditions reflect in the value of r [2sj ]. 
and can be such that this parameter is substantially 
smaller. Dynamics is then closer to that of the KS equa- 
tion. The intermediate coarsening regime narrows, and 
is followed by kinetic roughening. A surface morphology 
produced in these conditions [r — 0.5] is shown in Fig. 
\Hp), which can be compared with an AFM scan [(d)] of 
a Si target irradiated as in Jjjj. Again, agreement is ex- 
cellent. Note that the morphology now differs apprecia- 
bly from that of the KS equation, displayed in Fig.^c). 
While for Eq. (JSJ) a short-range ordered pattern coexists 
with long-range disorder and roughening, in the pure KS 
system disorder of the cellular structure is paradigmatic, 
see the height autocorrelations in Figs.^6), (c). Still, the 
time evolution of the roughness in Fig.|2f6) (o), predicted 
by Eq. (J5J for small r values, is similar to that of the KS 
case, Fig- Et&) (+): initial rapid growth is followed by 
much slower dynamics, and saturation to the stationary 
state. Such is also the experimental behavior found for 
nanostructures produced on Si, see Fig. 3 in 0(a). Com- 
paring the two plots in Fig. H[&), f° r small (non-zero) r 
values the small-scale nonlinearity A2 is seen to stabilize 
the linear instability earlier, and leads to smaller station- 
ary roughness. Moreover, in contrast with Fig.[2[a), Fig. 
Pfft) shows that for small or zero r values, the roughness 
does not have a local maximum as a function of time. 

In summary, we have introduced a continuum model 
for the formation of nanometric sized patterns by IBS. 
The model accounts within an unified framework for ex- 
perimental features of nanopatterns recently produced 
on diverse materials. Moreover, it leads to an effective 
interface equation providing new predictions. Thus, con- 
sidering dependencies Q on ion energy E of the features 
of the distribution of deposited energy, a, [i, a, the dot 
size l c behaves, in the large redeposition limit <j) < 1, as 
l c ~ [E + const.] 1 / 2 . For small E, this implies l c is energy 
independent, while l c ~ E 1 / 2 for large enough energies. 
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Observations exist |3 HH compatible with such energy 
dependence, although a systematic study assessing the 
importance of redeposition would be highly desirable. 

From a fundamental point of view, Eq. JSJ) also leads to 
new results. Specifically, this is a height equation with lo- 
cal interactions in which a pattern is stabilized with con- 
stant wavelength and amplitude, in contrast with conjec- 
tures for Id systems Although more theoretical work 
is still needed [e.g., regarding the asymptotic properties 
of Eq. (0)] this suggests that in 2d patterns, coarsening 
dynamics is indeed more complex than in Id |3l| . 
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FIG. 1: (a) Stationary state morphology from the numer- 
ical solution of Eq. JKJ with v = 2K, = 2, Ai = 0.1, 
A2 = 0.5, lateral size L = 256. Units are arbitrary. In- 
set: 2d-autocorrelation, showing high degree of short-range 
hexagonal order, (b) Same as before, for Ai = 1, L = 512. (c) 
same as (b) for the KS equation, A2 = 0. (d) 3x3 /j,m 2 AFM 
scan of a Si target irradiated as in Ref. [3 for 4 h. 
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FIG. 2: (a) W(t) and £ c (i) from Eq. © for Fig. [Ha). Dotted 
line in the linear (coarsening) region grows as an exponential 
(as t - 27 ). (b) W(t) for Fig.|D» (o), and Fig. Etc) (+). 



